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Abstract 

Seesaw-type and low-scale models of neutrino masses are reviewed, along 
with the corresponding structure of the lepton mixing matrix. The status of 
neutrino oscillation parameters as of June 2006 is given, including recent fluxes, 
as well as latest SNO, K2K and MINOS results. Some prospects for the next 
generation of experiments are given. This writeup updates the material presented 
in my lectures at the Corfu Summer Institute on Elementary Particle Physics in 
September 2005. 
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1 Introduction 

The historic discovery of neutrino oscillations [1-5] marks a turning point in particle 
and nuclear physics and implies that neutrinos have mass. This possibility has been 
first suggested by theory since the early eighties, both on general grounds and on the 
basis of different versions of the seesaw mechanism [6-11]. 

The general characterization of neutrino mass theories in SU(3) ® SU(2) ® U(l) 
terms provided a model-independent basis to analyse the seesaw [9,10]. It also indicated 
a fundamental difference between the lepton and the quark mixing matrices, namely, 
the appearance of new phases associated to the Majorana nature of neutrinos [9]. 

Irrespective of what the ultimate origin of neutrino mass may turn out to be, the 
basic gauge theoretic mechanism to account for the smallness of neutrino mass is in 
terms of the feebleness of B-L violation. The seesaw is one realization of the idea, by far 
not unique. There are two classes of theories of neutrino mass, that differ by the scale 
at which L-symmetry is broken. They are summarized in Sec. H3 The corresponding 
structure of the lepton mixing matrix that follows from theory is described Sec. EJ This 
forms the basis for the analysis of the data from current neutrino oscillation experi- 
ments [1-5] [12,13]. The status of neutrino mass and mixing parameters as determined 
from the world's neutrino oscillation data within the simplest CP-conserving three- 
neutrino mixing scheme is summarized in Sec. 15.11 [14] . In addition a determination of 
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the solar angle #i 2 , the atmospheric angle 6 2 3 and the corresponding mass squared split- 
tings Amg 0L and Am^ TM , one gets a constraint on the last angle in the three-neutrino 
leptonic mixing matrix, #13. Together with the small ratio Amg 0L /Am^ TM the angle #13 
holds the key for further progress in neutrino oscillation searches. 

Some attempts at predicting neutrino masses and mixing are given in Sec. 15.21 
Lepton number violating processes such as neutrinoless double beta decay [15, 16] are 
briefly discussed in Sec. 15.31 Searching for Ovflfl constitutes a very important goal 
for the future, as this will probe the fundamental nature of neutrinos, irrespective of 
the process that induces it, a statement known as the "black-box" theorem [17]. In 
addition, Qu/3(3 will be sensitive to the absolute scale of neutrino mass and to CP 
violation induced by the so-called Majorana phases [9], inaccessible in conventional 
oscillations [18-20]. Finally in Sec. El the robustness of the oscillation interpretation 
and the role of non-standard neutrino interactions in future precision oscillation studies 
is briefly mentioned. 



2 Dirac and Majorana masses 

Electrically charged fermions must be Dirac type. In contrast, electrically neutral 
fermions, like neutrinos (or supersymmetric " inos"), are expected to be Majorana- 
type on general grounds, irrespective of how they acquire their mass. Phenomenological 
differences between Dirac and Majorana neutrinos are tiny for most processes, such as 
neutrino oscillations: first because neutrinos are known to be light and, second, because 
the weak interaction is chiral, well described by the V-A form. 

The most basic spin 1/2 fermion corresponding to the lowest representation of the 
Lorentz group is given in terms of a 2-component spinor p, with the following free 
Lagrangean [9] 

C M = -ip^a^p - —p T a 2 p + H.C. (1) 

where <Ji are the usual Pauli matrices and 04 = —il, I being the 2x2 identity matrix. 
I use Pauli's metric conventions, where a. b = a^b^ = a ■ b + a 4 6 4 , a 4 = ia®. Under a 
Lorentz transformation, x — > Ax, the spinor p transforms as p —>■ S(A)p(A~ 1 x) where 
S obeys 

SV^S = A„ v a v (2) 

The kinetic term in Eq. is clearly invariant, and so is the mass term, as a result 
of unimodular property det 5 = 1. However, the mass term is not invariant under a 
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phase transformation 

P - e> (3) 
The equation of motion following from Eq. (JJJ is 

-i&n d nP = ma 2 p* (4) 

As a result of the conjugation and Clifford properties of the cr-matrices, one can verify 
that each component of the spinor p obeys the Klein-Gordon wave-equation. 

Start from the usual Lagrangean describing of a massive spin 1/2 Dirac fermion, 
given as 

C D = -#7/^ - m ( 5 ) 

where by convenience we use the chiral representation of the Dirac algebra 7/^7^+7^7^ = 
2 8^ v in which 75 is diagonal, 

In this representation the charge conjugation matrix C obeying 

C T = -C (7) 
C ] = C- 1 (8) 
C- 1 1, C = - 7 J (9) 

is simply given in terms of the basic conjugation matrix 02 as 

-(?:)■ 

In order to display clearly the relationship between the Majorana theory in Eq. ([TJ and 
the familiar Dirac Lagrangean in Eq. (J5)l. one splits a Dirac spinor as 

*D= [ ^ i ■ I I 

so that the corresponding charge-conjugate spinor ^/ C D = C ^ is the same as but 
exchanging and i- e. 

*!>=( ^ J • (12) 
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A 4-component spinor is said to be Majorana or self-conjugate if \I/ = C^ T which 
amounts to setting x — 0- Using Eq. we can rewrite Eq. (J5J) as 

2 2 
£d = -i Y pi a AP* Pl^pa + H.C. (13) 

a=l a=l 

where 

4> = ^{p2-ipx) (14) 

are the left handed components of ^> e> and of the charge-conjugate field ty c D , re- 
spectively. This way the Dirac fermion is shown to be equivalent to two Majorana 
fermions of equal mass. The U(l) symmetry of the theory described by Eq. (jSj) under 
^ d —* e ia ^D corresponds to continuous rotation symmetry between p\ and p2 

pi — > cos 6 pi + sin 9p2 
P2 — > — sin 9p\ + cos 6p2 

which result from the mass degeneracy between the p's, showing that, indeed, the 
concept of fermion number is not basic. 

The mass term in Eq. ((TJ vanishes unless p and p* are anti-commuting, so the 
Majorana fermion is, right from the start, a quantized field. The solutions to Eq. 
are easily obtained in terms of those of Eq. which are well known: 



r 2 

* M = (2vr)- 3 / 2 / d 3 kJ2(^) 1/2 [e lk - X Mk)u L r(k) +e- ik - x At(k)v Lr (k)}, 

J r =l 



(15) 



where u = C v T and E{k) = (k 2 + m 2 ) 1 ^ 2 is the mass-shell condition. The creation and 
annihilation operators obey canonical anti-commutation rules and, like the w's and v's, 
depend on the momentum k and helicity label r. The expression in Eq. (|T5|l describes 
the basic Fourier expansion of a massive Majorana fermion. It differs from the usual 
Fourier expansion for the Dirac spinor in Eq. (|16|) in two ways: 

• spinors are two-component, as there is a chiral projection on the w's and v 's 

• there is only one Fock space, particle and anti-particle coincide, showing that 
a massive Majorana fermion corresponds to one half of a conventional massive 
Dirac fermion. 
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The li's and u's are the same wave functions that appear in the Fourier decomposition 
the Dirac field 

r 2 

q D = (2tt)- 3 / 2 / d 3 k J2{^) 1/2 [e ik - x a r {k)u r {k) + e - lk - x bl{k)v r {k)] . (16) 

Using the helicity eigenstate wave-functions, 

ff. kuf(k) = ± | k | uf(k) (17) 
a- kvf(k) = =f | k \ vf(k) (18) 

one can show that, out of the 4 linearly independent wave functions u L (k) and v^(k), 
only two survive as the mass approaches zero, namely, u~[{k) and vj i {k) [21]. This 
way the Lee- Yang two-component massless neutrino theory is recovered as the massless 
limit of the Majorana theory. 

Two independent propagators follow from Eq. (JTJ), 

< | p(x) p*(y) | >= ia^Apix - y; m) (19) 
< | p(x) p(y) | >= m cr 2 Ap(x — y; m) (20) 

where Ap{x — y;m) is the usual Feynman function. The first one is the "normal" 
propagator that intervenes in total lepton number conserving (AL = 0) processes, 
while the one in Eq. (120)1 describes the virtual propagation of Majorana neutrinos in 
AL = 2 processes such as neutrinoless double-beta decay. 

The Lagrangean in Eq. (P) can easily be generalized to an arbitrary number of 
Majorana neutrinos, giving 
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L M = ~iJ2 pl a v d nP* M ^pl a 2pf3 + H.C. (21) 

a=l a,(3=l 

where the sum runs over the "neutrino-type" indices a and j3. By Fermi statistics the 
mass coefficients M a p must form a symmetric matrix, in general complex. This matrix 
can always be diagonalized by a complex nxn unitary matrix U (See [9] for the proof) 

M diag = U T MU . (22) 

When M is real its diagonalizing matrix U may be chosen to be orthogonal and, in 
general, the mass eigenvalues can have different signs. These may be assembled as a 
signature matrix 

77 = diag(+, +,..., -,-,..) (23) 
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S(7(3) ® 6(7(2) <g> (7(1) 


L a = {yai la) 


(1,2,-1) 




(1,1,2) 


/O (A \T 


(3,2,1/3) 


K 


(3, 1, -4/3) 


da 


(3,1,2/3) 


$ 


(1,2,1) 



Table 1: Matter and scalar multiplets of the Standard Model (SM) 



For two neutrino types there are two classes of models, one with rj = diag(+, — ) and 
another characterized by 77 = diag(+, +). The class with rj = diag(+, — ) contains as 
a limit the case where the two fermions make up a Dirac neutrino. Note that one can 
always make all masses positive by introducing appropriate phase factors in the wave 
functions, such as the factors of i in Eq. (|14jl . However, when interactions are added 
these signs become physical. As emphasized by Wolfenstein, they play an important 
role in the discussion of Qv(3(3 (neutrinoless double beta decay) [22]. 



3 The origin of neutrino mass 

Table ^ gives the fifteen basic building blocks of matter. They are all 2-component 
sequential "left-handed" chiral fermions, one set for each generation. Parity violation in 
the weak interaction is incorporated "effectively" by having "left" and "right" fermions 
transform differently with respect to the SU (%)®SU(2)®U{1) gauge group. In contrast 
to charged fermions, neutrinos come only in one chiral species. 

It has been long noted by Weinberg [8] that one can add to the Standard SU{2>) ® 
SU{2) ® £7(1) Model (SM) an effective dimension-five operator = AL$L$ where L 
denotes a lepton doublet for each generation and $ is the SM scalar doublet. 




Figure 1: Dimension five operator responsible for neutrino mass [8]. 
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Once the electroweak symmetry breaks through the nonzero vacuum expectation 
value (vev) ($), Majorana neutrino masses oc ($) 2 are induced, in contrast to the 
masses of the charged fermions which arise from basic renormalizable interactions, and 
are linear in ($). Moreover, the dimension-five operator violates lepton number by 
two units (AL = 2), whereas the charged fermion masses arise from renormalizable L- 
conserving Yukawa interactions. This naturally accounts for the smallness of neutrino 
masses irrespective of the specific origin of neutrino mass. From such general point of 
view the emergence of Dirac neutrinos would be a surprise, justified only in the presence 
of an "accidental" lepton number symmetry. For example, neutrinos could naturally 
get very small Dirac masses via mixing with a bulk fermion in models involving extra 
dimensions [23-25]. Barring such very special circumstances, gauge theories give rise 
to Majorana neutrinos. 

Little more can be said from first principles about the mechanism giving rise to the 
operator in Fig. its associated mass scale or its Savour structure. For example, the 
strength A of the operator may be suppressed by a large scale M x in the denominator 
(top-down) scenario, leading to 



where Ao is some unknown dimensionless constant. Gravity, which in a sense "belongs" 
to the SM, could induce the dimension-five operator 0, providing the first example of 
a top-down scenario with Mx = Mp, the Planck scale [26]. In this case the magnitude 
of the resulting Majorana neutrino masses are too small. 

Alternatively, the strength A of the operator may be suppressed by small param- 
eters (e.g. scales, Yukawa couplings) in the numerator and/or loop-factors (bottom-up 
scenario). Both classes of scenarios are viable and have many natural realizations. 
While models of the top-down type are closer to the idea of unification, bottom-up 
schemes are closer to experimental test. 

Models of neutrino mass may also be classified according to whether or not addi- 
tional neutral heavy states are present, in addition to the three isodoublet neutrinos. 
As an example, such leptons could be SU(3) <E> SU(2) <8> U(l) singlet "right-handed" 
neutrinos. In what follows we first consider top-down, then bottom-up scenarios. 

3.1 Seesaw and related models 

The most popular top-down scenario is the seesaw [6] . The idea is to generate the oper- 
ator by the exchange of heavy states, both fermions (type-I) and scalars (type-II), as 
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shown in Fig. |21 This can be implemented in many ways, with different gauge groups 
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Figure 2: Two types of seesaw mechanism [6-11] 



and multiplet contents. The main point is that, as the masses of the intermediate states 
go to infinity, neutrinos naturally become light [7]. The seesaw provides a simple re- 
alization of Weinberg's operator [8]. Note that the seesaw idea does not require the 
gauging of B-L, nor does it require it to be broken spontaneously. In fact, most of the 
physics it encodes, and which has been brilliantly confirmed by the recent oscillation 
experiments, lies in its "effective" low-energy form [9]. Only in low-scale schemes like 
the inverse-seesaw discussed in Sec. IH.2.HI effects associated with "seesaw-dynamics" 
may be observable (see Sees. 15.4.11 and 15.4. 2j) . 

3.1.1 "Effective" seesaw 

Much of the low energy phenomenology, such as that of neutrino oscillations is blind 
to the details of the underlying seesaw theory at high energies, e. g. its gauge group, 
multiplet content or the nature of B-L. For this purpose the most general way to 
describe the physics of the seesaw is to characterize it, effectively, in terms of the SM 
gauge structure [9]. In the basis vl, i/£, corresponding to the three "left" and three 
"right" neutrinos, respectively, the seesaw mass matrix has SU(2) triplet, doublet and 
singlet terms described as [9] 



Here we use the original notation of reference [9], where the "Dirac" entry is propor- 
tional to (<&), the Mi comes from a triplet vev, and M2 may be added by hand, as it 
is a gauge singlet. The particular case Mi = was first mentioned in Ref. [6]. 

Note that, though symmetric, by the Pauli principle, the matrix Ai v is complex, so 
that its Yukawa coupling sub-matrices D as well as Mi and M 2 are also complex matri- 
ces, the last two symmetric. It is diagonalized by performing a unitary transformation 




(24) 
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u v , 

6 

Vi = y](U v )ian a , (25) 

a=l 

so that 

U?MvU v = tiag{m i ,M i ). (26) 

This yields 6 mass eigenstates, including the three light neutrinos with masses m^, and 
three two-component heavy leptons of masses Mj. The light neutrino mass states i/j 
are given in terms of the flavour eigenstates via eq. (|25j) . The effective light neutrino 
mass, obtained this way is of the form 

m„ ~ Mi - DM 2 ~ 1 D T . (27) 

The smallness of light neutrino masses is understood by assuming M 2 3> D ^> Mi. 
The above general structure forms the basis for the description of the seesaw lepton 
mixing matrix [9], given in Sec. 14.31 

While it constitutes the most general description, and also the common denominator 
of all seesaw schemes, such an "effective" seesaw does not give a dynamical insight on 
the origin of neutrino mass. For this reason we now turn to schemes where lepton 
number symmetry is broken spontaneously. 



3.1.2 The "1-2-3" seesaw mechanism 

The simplest possibility for the seesaw is to have ungauged lepton number. It is also the 
most general, as it can be studied in the framework of just the SU(3) <g> SU{2) ® U(l) 
gauge group. The mass terms in eq. (J24f are given by triplet, doublet and singlet vevs, 
respectively, as [10] 

Y 3 v 3 Y v ($) 
Y/ ($) Y lVl 

As already mentioned, the Yukawa coupling sub-matrices Y v as well as Y 3 and Y 1 are 
complex matrices, the last two symmetric. 

The new dynamical insight provided by such "1-2-3" seesaw containing singlet, 
doublet and triplet scalar multiplets, is that they obey a simple vev seesaw relation of 
the type 

v 3 Vi ~ v 2 2 with v i ^> v 2 ^> t>3 (29) 

where v 2 = ($) denotes the SM Higgs doublet vev, fixed by the W-boson mass. This 
hierarchy implies that the triplet vev v 3 — >• as the singlet vev v% grows. This is 
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consistent with the minimization of the corresponding SU (3) ® SU(2) ® U(l) invariant 
scalar potential, and implies that the triplet vev is "induced" . 

Small neutrino masses arise either by heavy SU(3) ® SU(2) ® U(l) singlet "right- 
handed" neutrino exchange (type I) or by the small effective triplet vev (type II), as 
illustrated in Fig. 121 The effective light neutrino mass becomes, 



The corresponding seesaw diagonalization matrices can be given explicitly as a system- 
atic matrix perturbation series expansion in DM^ 1 , given in Ref. [10]. 

In such "1-2-3" seesaw, since lepton number is ungauged, there is a physical Gold- 
stone boson associated with its spontaneous breakdown, the majoron [27]. Its profile 
can be determined on symmetry grounds and its couplings to neutrinos can be found 
systematically, see [10] for details. 

3.1.3 Left-right symmetric and SO (10) 

A more symmetric setting for the seesaw is a gauge theory containing B-L as a gen- 
erator, such as SU(3) ® SU{2) L ® SU{2) R <g> U{\)b-l or the unified models based on 
50(10) or £"(6) [6,7,11]. For example in 50(10) each matter generation is naturally 
assigned to a 16 (spinorial in 5*0(10)) so that the 16 . 16 . 10 and 16 . 126 . 16 
couplings generate all entries of the seesaw neutrino mass matrix, 



Here the basis is v L , i/£, as before, Y L and Y R denote the Yukawas of the 126 of 5*0(10), 
whose vevs (Al,r) give rise to the Majorana terms. They correspond to Y\ and Y3 of 
the simplest "1-2-3" model. On the other hand Y v denotes the 16 . 16 . 10 Dirac 
Yukawa coupling. In 5O(10) one has a discrete parity symmetry which implies Yl = Yr 
and Y v = Yj as recently emphasized in Ref. [28]. Since this may get broken, we prefer 
to keep, for generality, Yl, Yr as independent. 

Small neutrino masses are induced either by heavy SU(3) <g> SU{2) ® U{1) singlet 
"right-handed" neutrino exchange (type I) or heavy scalar boson exchange (type II) as 
illustrated in Fig. El The matrix M. v is diagonalized by a unitary mixing matrix U v 
as before. The diagonalization matrices can be worked out explicitly as a perturbation 
series, using the same method of Ref. [10]. This means that the explicit formulas for 
the 6x6 unitary diagonalizing matrix U given in Ref. [10] also hold in the left-right 
case, provided one takes into account that v% — > (Ar) and v% — > (A^). 



m v ~ Y 3 v 3 - Y V Y X - V Y V 



T 



<$> 2 



(30) 




(31) 



11 



The effective light neutrino mass, obtained this way is of the form 



m v « Y L (A L ) - Y V Y R Y v 



t(*) 2 



(32) 



(A*) 



We have the new vev seesaw relation 



(A L ) (A fl ) ~ ($) 2 



(33) 



which naturally follows from minimization of the left-right symmetric scalar potential, 
together with the vev hierarchy 



This implies, as before, that both type I and type II contributions vanish as (Ar) — > oo. 
Notice that, strictly speaking, this version of left-symmetry is inconsistent with type-I 
seesaw [7] and requires the full form of the seesaw mass matrix [9,10]. There are, 
however, tripletless left-right seesaw variants where a similar vev-seesaw formula holds, 



If one can arrange for the breakdown of parity invariance to be spontaneous, then 
the smallness of neutrino masses gets correlated to the observed maximality of parity 
violation in low-energy weak interactions, as stressed by Mohapatra and Senjanovic [7]. 
However elegant this connection may be, it is phenomenologically not relevant, in view 
of the large value of the B-L scale. The latter is required both to fit the neutrino masses, 
as well as to unify the gauge couplings. Another important difference with respect to 
the simplest "1-2-3" seesaw is the absence of the majoron, which is now absorbed as 
the longitudinal mode of the gauge boson corresponding to the B-L generator which 
picks up a huge mass. 

3.1.4 "Double" seesaw mechanism 

One can add any number of (anomaly-free) gauge singlet leptons Si to the SM, or any 
other gauge theory [9]. For example, in 50(10) and E(6) one may add leptons outside 
the 16 or the 27, respectively. New important features may emerge when the seesaw 
is realized with non-minimal lepton content. Here we mention the seesaw scheme 
suggested in Ref. [29] motivated by string theories [30]. The model extends minimally 
the particle content of the SM by the sequential addition of a pair of two-component 
SU(3) ® SU(2) ® 17(1) singlet leptons, z>f, Si, with i a generation index running over 



(A L ) « ($) « (A R ) 



(34) 



see Refs. [31-34] and Sec. l3~T31 
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1, 2, 3. In the v, z/ c , S basis, the 9x9 neutral leptons mass matrix M. v is given as 

/ YJ ($) 

Yv ($> 
M 



My 



\ 



(35) 



\ 

M T , 

/* / 

in the basis vl, v c Ll Sl. Again Y v is an arbitrary 3x3 complex Yukawa matrix, M 
and \i are SU(2) singlet complex mass matrices, \x being symmetric. Notice that it has 
zeros in the Pl~ v l and v c L -v c L entries, a feature of several string models [30]. 

For fi ^> M one has to first approximation that the Si decouple leaving the simpler 
seesaw at scales below that. In such a "double" seesaw scheme the three light neutrino 
masses are determined from 

m v = M T_I fxM- x Y u , (36) 

The mass generation is illustrated in Fig. |3J A new feature is that there are two 




V 



Figure 3: "Double" and "inverse" seesaw mechanism [29]. 

independent scales, \x and M, of which only \x breaks the B-L symmetry. Both scales 
can be large but, as we will see in Sec. 13.2.31 it is natural for fi to be small [29], instead 
of large, /i <C M, see Sec. 13.2.31 For the case /x = M one formally recovers the usual 
seesaw form, this is useful to present some results in simplified form, as used in Fig. ITU 
Irrespective of what sets its scale, the entry /i may be proportional to the vev of 
an SU(3) ® SU(2) ® U{\) singlet scalar, in which case the model contains a singlet 
majoron, see Sec. 15.4.11 

3.1.5 "Novel" SO (10) seesaw mechanism 

The seesaw is a mechanism which allows for many possible realizations. Schemes 
leading to the same pattern of neutrino masses may differ in many other respects. 
Correspondingly, there are many types of seesaw. In addition to type I [6,7] [11], type 
II [9,10], there are extended seesaw schemes, like type-Ill [31-33] and the double/inverse 
seesaw described above. 
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Here I turn to yet another seesaw that has recently been suggested in Ref. [34]. 
It belongs to the class of supersymmetric SO(10) models with broken D-parity. In 
addition to the states in the 16, it contains three sequential gauge singlets Sn with 
the following mass matrix 



/ F( XL )\ 

Yj ($) F ( X r) 
\F t ( Xl ) F t ( Xr ) / 



(37) 



in the same basis u L , z/£, Sl as previously. Notice that it has zeros along the diagonal, 
specially in the vl-vl and v L -v1 entries, thanks to the fact that there is no 126, a 
feature of several string models [29,30]. The resulting light neutrino mass matrix is 

<2 



M 



unif 



Y V (FF- V ) T + {FF- L )Y t 



T 



(3? 



where M un if is the unification scale, F and F denote independent combinations of 
Yukawa couplings of the Sn- One can see that the neutrino mass is suppressed by 
the unification scale M un ;f independently of the B-L breaking scale. In contrast to all 
previous seesaws, this one is linear in the Dirac Yukawa couplings Y u , as illustrated 
in Fig. |3J It is rather remarkable that one can indeed take the B-L scale as low as 



$ \ / XR 



i 

\XL 
I 

-* 



I 



Xr 



: $ 



s 



Figure 4: Low B-L scale 50(10) seesaw mechanism [34]. 



TeV without generating inconsistencies, neither with neutrino masses, nor with gauge 
coupling unification [34]. 



3.2 Low-scale models 

There are many models of neutrino mass where the operator is induced from physics 
at accessible scales, TeV or less. The smallness of its strength may be naturally achieved 
due to loop and Yukawa couplings suppression. Moreover, the strength of the operator 
may be suppressed by small lepton number violating parameters that appear in its 
numerator, instead of its denominator, as commonly assumed. The latter correspond 



14 



to the seesaw-type schemes previously discussed. The former correspond to the class 
we are about to describe. Before I do so, let me emphasize that such models are also 
"natural" in t'Hooft's sense [35]: "an otherwise arbitrary parameter may be taken as 
small when the symmetry of the Lagrangean increases by having it vanish". 

Since all particles are at the TeV scale, these models naturally lead to possibly 
testable phenomenological implications, including lepton flavour violation and modifi- 
cations in muon and tau decays. 

Moreover, when the breaking of lepton number entailed in these models takes place 
spontaneously, the corresponding Goldstone boson has a remarkable property: it may 
couple substantially to the SM Higgs boson, which can therefore have a sizeable invisible 
decay branching ratio [36-39] 



where J is the majoron. This show that, although neutrino masses are small, the 
neutrino mass generation may have very important implications for the mechanism of 
electroweak symmetry breaking. One must therefore take into account the existence 
of the invisible channel in designing Higgs boson search strategies at future collider 
experiments [40,41]. Further discussion in Sec. 15.4. Tl 

3.2.1 Radiative models 

Neutrino masses may be induced by calculable radiative corrections [42]. For example, 
they can arise at the two-loop level [43] as illustrated in Fig. |3J Up to a logarithmic 
factor one has, schematically, 



in the limit where the doubly-charged scalar k is much heavier than the singly charged 
one. Here I denotes a charged lepton, / and h are their Yukawa coupling matrices and 
Yi denotes the SM Higgs Yukawa couplings to charged leptons. Here (a) denotes an 
SU (3) ® SU (2) ® Z7 (1) singlet vev used in Ref. [44] . Clearly, even if the proportionality 
factor Ao is large, the neutrino mass can be made naturally small by the presence of a 
product of five small Yukawas and the appearance of the two-loop factor. 

3.2.2 Supersymmetry and neutrino mass 

Low energy supersymmetry can be the origin of neutrino mass [45]. The intrinsically 
supersymmetric way to break lepton number is to break the so-called R parity. This 
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Figure 5: Two- loop origin for neutrino mass [43,44]. 



could happen spontaneously, driven by a nonzero vev of an SU (3)®SU (2)®Z7(1) singlet 
sneutrino [46-48], and leads to an effective model characterized by purely bilinear R 
parity violation [49] . This also serves as reference model, as it provides the minimal way 
to add neutrino masses to the MSSM, we call it RMSSM. Neutrino mass generation 
takes place in a hybrid scenario, with one scale generated at tree level and the other 
induced by "calculable" radiative corrections [50]. Here the two blobs in the graph 



Fig. El denote AL = 1 insertions, while the crossed blob accounts for chirality flipping, 
and A denotes the trilinear soft supersymmetry breaking coupling. The general form 
of the expression is quite involved but the approximation 



holds in some regions of parameters. The neutrino mass spectrum naturally follows a 
normal hierarchy, with the atmospheric scale generated at the tree level and the solar 
mass scale arising from calculable loops. 

3.2.3 "Inverse" seesaw mechanism 

Here we mention that there are also low-scale tree-level neutrino mass schemes with 
naturally light neutrinos. One is the inverse seesaw scheme suggested in [29] and 
already described in Sec. 13.1.41 The mass matrix is the same as that of the double 




Figure 6: Loop origin of solar mass scale in RMSSM [50]. 




(41) 
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seesaw model, given in eq. (|35|) . The only difference is that now the entry \i is taken 
very small, e. g. ll <C Y v (0) <C M [29]. Notice that for small ll neutrino masses vanish 
with li, 

m u = (<S>} 2 Yj M T_1 iiM- x Y v 

as illustrated in Fig. 01 

The fact that the neutrino mass vanishes as ll — > is just the opposite of the 
behaviour of the seesaw formulas in Eqs. (J2*7|) (jHOJ) and (JH2J); thus this is sometimes 
called inverse seesaw model of neutrino masses. The entry ll may be proportional 
to the vev of an SU(2) singlet scalar, in which case spontaneous B-L violation leads 
to the existence of a majoron [51]. This would be characterized by a relatively low 
scale, so that the corresponding phase transition could take place after the electroweak 
transition. The naturalness of the model stems from the fact that in the limit when 
Li — > lepton number is recovered, increasing the symmetry of the theory. 

4 The lepton mixing matrix 

In any gauge theory in order to identify physical particles one must diagonalize all 
relevant mass matrices, which typically result from gauge symmetry breaking. Mech- 
anisms giving mass to neutrinos generally imply the need for new interactions whose 
Yukawas (like Y v in seesaw-type schemes) will coexist with that of the charged leptons, 
Y\. Since in general these are independent one has that, like quarks, massive neutri- 
nos generally mix. The structure of this mixing is not generally predicted from first 
principles. Whatever the ultimate high energy gauge theory may be it must be broken 
to the SM at low scales, so one should characterize the structure of the lepton mixing 
matrix in terms of the SU(3) <8> SU (2) ® U (1) structure. The procedure is the familiar 
one from the quark sector. 

4.1 Dirac case 

Here we derive the structure of the lepton mixing matrix of massive Dirac neutrinos 
and its parametrization, as presented in Ref. [9]. From the start V can be parametrized 

as 

n 
i<j 
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where 



w„(7) = expi(^7oA*) 



0=1 



is a diagonal unitary "Cartan" matrix, described by n — 1 real parameters 7 a On 
the other hand each factor 



^ab 



(r) ab ) = exp Ys^abAl - r]* ab A a b ) 



0=1 



is a complex rotation in ab with parameter rj ab = |?y a b|exp 2# a ft. For example, 

/ cia e^sia ... \ 



W12 (?7 



12 



2 Sl2 C12 



V 



(43) 



/ 



Once the charged leptons and Dirac neutrino mass matrices are diagonal, one can still 
rephase the corresponding fields by u (a) and ajQ^ — a), respectively, keeping invariant 
the form of the free Lagrangean. This results in the form 



V = u (a) Yluijfaj) ul 



0K a ) , 



(44) 



i<j 



where we are still free to choose the n — 1 a-values associated to Dirac neutrino phase 
redefinitions. Using the conjugation property 



a;o(a)^ a6 (|^ ai ,|exp iO ab ) a;J(a) = u; afl [|?7 ab |exp i(a a + 6 ab - a b )) 



(45) 



we arrive at the final Dirac lepton mixing matrix which is, of course, identical in form 
to that describing quark mixing. It involves a set of 



n[n 



l)/2 mixing angles 0^ and n(n — l)/2 — (n — 1) CP phases. (46) 



where (n — 1) phases were eliminated. This is the parametrization as originally given in 

Ref. [9], with unspecified factor ordering. From Eq. fJ4fij) one sees that for n = 3 there 

are 3 angles and precisely one leptonic CP violating phase, just as in the Kobayashi- 

Maskawa matrix describing quark mixing. Two of the three angles are involved in solar 

choosing an overall relative phase between charged leptons and Dirac neutrinos we can take 
V as unimodular, i. e. det V = 1, so that J2a=i 7a = 1- 
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and atmospheric oscillations, so we set # 12 = # SOL and #23 = #atm- The last angle in the 
three-neutrino leptonic mixing matrix is #13, 

/ da e^s 13 \ 
u 13 = 10 . (47) 

V -e-^s 13 da / 

A convenient ordering prescription is to take 23 x 13 x 12, or "atmospheric" x 
"reactor" x "solar", with the phase being associated to #13. In summary, if neutrinos 
masses are added a la Dirac their charged current weak interaction has exactly the 
same structure as that of quarks. 

4.2 Majorana case 

Here we consider the form of the lepton mixing matrix in models where neutrino masses 
arise in the absence of right-handed neutrinos, such as those in Sec. IH.2.11 The unitary 
form also holds, to a good approximation, to models where SU(2) doublet neutrinos 
mix only slightly with other states, like high-scale seesaw models. 

For n generations of Majorana neutrinos the lepton mixing matrix has exactly the 
same form given in Eq. f!44j) . The difference is that in the case of Majorana neutrinos 
their mass term is manifestly not invariant under rephasings of the neutrino fields. As 
a result, the parameters a in Eq. f!44j) can not be used to eliminate n — 1 Majorana 
phases as we just did in Sec. 14.11 Consequently these are additional phases [9] which 
show up in L-violating processes [19,20]. Such new sources of CP violation in gauge 
theories with Majorana neutrinos are called "Majorana phases". They already exist 
in a theory with just two generations of Majorana neutrinos [9], n = 2, whose mixing 
matrix is described by 

\ -e l9l2 s u d 2 / 

where 0i2 is the Majorana phase (recall that Cabibbo mixing has no CP phase). Such 
"Majorana" CP phases are, in a sense, mathematically more "fundamental" than the 
Dirac phase, whose existence, as we just saw, requires at least three generations. 

For the case of three neutrinos the lepton mixing matrix can be parametrized as [9] 

K = ^23^13^12 (49) 

where each factor in the product of the cj's is effectively 2x2, characterized by an angle 
and a CP phase. Such symmetrical parameterization of the lepton mixing matrix, K 
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can be written as: 



c 12 c 13 s 12 c 13 e^ s 13 e 1 ^ 

-si 2 c 23 e-^ 12 - c^s^se^-^ c 12 c 23 - s 12 s 13 s 23 e i ^ ia+ * B -^ Ci 3 s 23 e ife 
. si2S 23 e-^ 12+te) - c 12 s 13 c 23 e-^« -c 12 s 23 e-^ _ Sl2Sl3 c 23e ! (*^fe) Cl3C23 

All three CP phases are physical [19]: 012,023 and 0i 3 . The "invariant" combination 
5 = 0i2 + </> 23 — 0i3 corresponds to the "Dirac phase" . If neutrinos are of Dirac type, 
only a single phase (say 0i 3 ) may be taken to be non-zero. This phase corresponds to 
the phase present in the Kobayashi-Maskawa matrix, and this is the one that affects 
neutrino oscillations. The other two phases are associated to the Majorana nature 
of neutrinos and show up only in lepton-number violating processes, like neutrinoless 
double beta decay [19,20]. 

An important subtlety arises regarding the conditions for CP conservation in gauge 
theories of massive Majorana neutrinos. Unlike the case of Dirac fermions, where CP 
invariance implies that the mixing matrix should be real, in the Majorana case the 
condition is [21] 

K* = Krj 

where rj = diag(+, +, , , ..) is the signature matrix describing the relative signs of the 
neutrino mass eigenvalues that follow from diagonalizing the relevant Majorana mass 
matrix, if one chooses to use real diagonalizing matrices [22]. Consequently say, for 
n = 2, both 0i2 = 7r/2 and i2 = correspond to CP conservation, as emphasized by 
Wolfenstein. These important signs determine the CP properties of the neutrinos and 
play a crucial role in Qvflfi . 

4.3 Seesaw-type mixing 

The most general theory of neutrino mass is effectively described in SM terms by (n, m), 
n being the number of ££7(3) <g> £77(2) ® £7(1) isodoublets and m/fl the number of 
SU(3) <S> ££7(2) ® £7(1) isosinglet two-component leptons (the case m = was given 
above). Here we assume an arbitrary number of gauge singlets, since they carry no 
anomaly. The usual seesaw in Sees. 13.1.11 13.1.21 and 13.1.31 has m = n = 3, while the 
extended seesaw in Sees. I3.1.4| 13.1.51 and 13.2.31 have m = 2n = 6. Isosinglets have in 
general a gauge and Lorentz invariant Majorana mass term. The procedure holds in 
any scheme of Majorana neutrino masses where isosinglet and isodoublet mass terms 
coexist [9]. 
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The effective form of such a "seesaw" lepton mixing matrix has, in addition to 
Majorana phases, many doublet-singlet mixing parameters, in general complex [9]. Its 
general structure is substantially more complex than "usual", being described by a 
rectangular matrix, called K. As a result one finds that leptonic mixing as well as CP 
violation may take place even in the massless neutrino limit [52,53]. 

The existence of these neutral heavy leptons could be inferred from low energy weak 
decay processes, where the neutrinos that can be kinematically produced are only the 
light ones. The mixing matrix describing the charged weak interactions of the light 
(mass-eigenstate) neutrinos is effectively non-unitary, since the coupling of a given light 
neutrino to the corresponding charged lepton is decreased by a certain factor associated 
with the heavy neutrino coupling. There are constraints on the strength of the such 
mixing matrix elements that follow from weak universality and low energy weak decay 
measurements, as well as from LEP. 

The full weak charged current mixing matrix K of the general (n, m) models involves 

n(n + 2m — l)/2 mixing angles 0y and n(n + 2m — l)/2 CP phases 0y . (50) 

For the explicit parametrization the reader is referred to the original paper, Ref. [9]. 
One sees that, for example, the usual seesaw model [labeled (3, 3) in our language] 
is characterized by 12 mixing angles and 12 CP phases (both Dirac and Majorana- 
type) [9]. 

This number far exceeds the corresponding number of parameters describing the 
charged current weak interaction of quarks. As already mentioned, the reason is 
twofold: (i) neutrinos are Majorana particles, their mass terms are not invariant under 
rephasings, and (ii) the isodoublet neutrinos in general mix with the 577(3) <g> SU(2) <g> 
£7(1) singlets. As a result, there are far more physical CP phases that may play a role 
in neutrino oscillations and/or leptogenesis (see below). 

Another important feature which arises in any theory based on SU(3) <g> £77(2) <g> 
£7(1) where isosinglet and isodoublet lepton mass terms coexist is that the leptonic 
neutral current is non-trivial [9]: there are non-diagonal couplings of the Z to the 
mass-eigenstate neutrinos. They are expressed as a projective hermitian matrix 

P = K ] K. 

This contrasts with the neutral current couplings of mass-eigenstate neutrinos in schemes 
where lepton number is conserved fSec. 14. lj) or where no isosinglet leptons are present, 
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i.e., m = (Sec. 14. 2|) . In both cases, just as for quarks, the neutral current couplings 
are diagonal (Glashow-Iliopoulos-Maiani mechanism). 

Before we close, note that, in a scheme with m < n, n — m neutrinos will remain 
massless, while 2m neutrinos will acquire Majorana masses, m light and m heavy [9]. 
For example, in a model with n = 3 and m = 1 one has one light and one heavy 
Majorana neutrino, in addition to the two massless ones. In this case clearly there will 
be less parameters than present in a model with m = n. 

5 Phenomenology 

Obviously the first phenomenological implication of neutrino mass models is the phe- 
nomenon of neutrino oscillations, required to account for the current solar and atmo- 
spheric neutrino data. The interpretation of the data relies on good calculations of the 
corresponding fluxes [54,55], neutrino cross sections and response functions, as well as 
on an accurate description of neutrino propagation in the Sun and the Earth, taking 
into account matter effects [56,57]. 

5.1 Status of neutrino oscillations 

Current neutrino oscillation data have no sensitivity to CP violation. Thus we neglect 
all phases in the analysis and take, moreover, the simplest unitary 3-dimensional form of 
the lepton mixing matrix in Eq. (I49j) with the three phases set to zero. In this approx- 
imation oscillations depend on the three mixing parameters sin 2 # 12 , sin 2 # 2 3> sin 2 # 13 
and on the two mass-squared splittings Am 2 OL = Am^ = m| — m 2 and Am 2 TM = 
Amgj = m\ — m\ characterizing solar and atmospheric neutrinos. The hierarchy 
^ m soL ^ ^ m ATM implies that one can set Am 2 OL = 0, to a good approximation, 
in the analysis of atmospheric and accelerator data. Similarly, one can set Am 2 TM to 
infinity in the analysis of solar and reactor data. 

The world's neutrino oscillation data and their analysis, as of June 2006, are given 
in Ref. [14] and will not be repeated here. The new developments are: new Standard 
Solar Model [58], new SNO salt data [59], latest K2K [60] and MINOS [61] data. 
These are briefly described in Appendix C of hep-ph/0405172 (v5). In what follows we 
summarize the updated results of the analysis which takes into account all these new 
data. Apart from the "positive" data already mentioned, the analysis also includes the 
constraints from "negative" oscillation searches at reactor experiments, CHOOZ and 
Palo Verde. 
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Figure 7: Three-neutrino regions allowed by the world's neutrino oscillation data at 
90%, 95%, 99%, and 3a C.L. for 2 d.o.f. as of June 2006, from Ref. [14]. In top panels 
Ax 2 is minimized with respect to undisplayed parameters. 

The three-neutrino oscillation parameters that follow from the global oscillation 
analysis in Ref. [14] are summarized in Fig. In the upper panels of the figure the 
Ax 2 is shown as a function of the three mixing parameters sin 2 # 12 , sin 2 8 2 3, sin 2 # 13 and 
two mass squared splittings Am\ x , Am^ , minimized with respect to the undisplayed 
parameters. The lower panels show two-dimensional projections of the allowed regions 
in the five- dimensional parameter space. In addition to a confirmation of oscillations 
with Am 2 TM , accelerator neutrinos provide a better determination of Am 2 TM . For 
example, comparing dashed and solid lines in Fig. one sees that the inclusion of the 
new data (mainly MINOS [61]) leads to a slight increase in Am 2 TM and an improvement 
on its determination (see [14] for details). On the other hand reactors [3] have played a 
crucial role in selecting large-mixing-angle (LMA) oscillations [62] out of the previous 
"zoo" of solutions [63]. The best fit values and the allowed 3a ranges of the oscillation 
parameters from the global data are summarized in Table [21 

Note that in a three-neutrino scheme CP violation disappears when two neutrinos 
become degenerate or when one of the angles vanishes [64]. As a result CP violation 
is doubly suppressed, first by a = Am 2 OL /Am 2 TM and also by the small mixing angle 
#13. The left panel in Fig. |H1 gives the parameter a, as determined from the global % 2 
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parameter 


Dest nt 


3c range 


Am 2 21 [10" 5 


eV 2 ] 


7.9 


7.1-8.9 


Am^ [10~ 3 


eV 2 ] 


2.6 


2.0-3.2 


sin 6 l i2 




0.30 


0.24-0.40 


sin 2 6> 23 




0.50 


0.34-0.68 


sin 2 6»i3 




0.00 


< 0.040 



Table 2: Neutrino oscillation parameters as of June 2006, from Ref. [14]. 




Figure 8: Determination of a = Am 2 OL /Am 2 TM and bound on sin 2 #13 from data, as of 
June 2006, from Ref. [14]. 

analysis. The right panel shows the impact of different data samples on constraining 
6*13. One sees that, although for larger Am 2 TM values the bound on sin 2 #13 is dominated 
by CHOOZ, this bound deteriorates quickly as Am 2 TM decreases (see Fig. IS]), so that 
the solar and KamLAND data become relevant. 

There is now a strong ongoing effort aimed at probing # 13 and CP violation in future 
neutrino oscillation searches at reactors and accelerators [65-67]. As we saw, the basic 
parameters a = Am 2 OL /Am 2 TM and sin 2 #13 characterizing the strength of CP violation 
in neutrino oscillations are small. Prospects for probing sin 2 # 13 at long baseline reactor 
and accelerator neutrino oscillation experiments are given in Ref. [68], 

Information on sin 2 # 13 may also come from a totally different class of studies of the 
day /night effect in large water Cerenkov solar neutrino experiments such as UNO or 
Hyper-K [69] [70]. 
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5.2 Predicting neutrino masses and mixing 

Gauge symmetry alone is not sufficient to predict particle mixings, neither for the 
quarks, nor for the leptons: such "flavour problem" has remained with us for a while. 

As we saw in Sec. 15. II five of the basic parameters of the lepton sector are currently 
probed in neutrino oscillation studies. These point towards a well defined pattern of 
neutrino mixing angles, quite distinct from that of quarks. Such pattern is not easy to 
account for in the context of unified schemes where quarks and leptons are connected. 
The data seem to indicate an intriguing complementarity between quark and lepton 
mixing angles [71-74]. 

There has been a rush of papers attempting to understand the values of the leptonic 
mixing angles from underlying symmetries at a fundamental level. For example the 
following form of the neutrino mixing angles has been proposed [75] 



tan $atm 


= tan 2 9° 23 


= 1 


sin 2 #chooz 


= sm 2 9° 13 


= 


tan 2 6> S0L 


= tan 2 6° 12 


= 0.5 



Such Harrison-Perkins-Scott pattern [76] could result from some kind of flavour symme- 
try, valid at a very high energy scale where the dimension-five neutrino mass operator 
arises. 

One approach to predict neutrino masses and mixing angles was the idea that 
neutrino masses arise from a common seed at some "neutrino mass unification" scale 
Mx [77] , very similar the merging of the SM gauge coupling constants at high energies 
due to supersymmetry [78]. However, in its simplest form this very simple theoretical 
ansatz is now inconsistent (at least if CP is conserved) with the current observed value 
of the solar mixing angle 9\2 inferred from current data. 

A more satisfactory and fully viable alternative realization of the "neutrino mass 
unification" idea employs an A\ flavour symmetry introduced by Ernest Ma, in the 
context of a seesaw scheme [79]. Starting from three-fold degeneracy of the neutrino 
masses at a high energy scale, a viable low energy neutrino mass matrix can indeed 
be obtained in agreement with neutrino data as well as constraints on lepton flavour 
violation in /i and r decays. The model predicts maximal atmospheric angle and 
vanishing # 13 , 

#23 = tt/4 and 9 13 = . 
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Although the solar angle 6*12 is unpredicted, one expects 2 

012 = 0(1)- 

When CP is violated #13 becomes arbitrary and the Dirac phase is maximal [81]. 

Within such flavour symmetric seesaw scheme one can show that the lepton and 
slepton mixings are intimately related. The resulting slepton spectrum must necessarily 
include at least one mass eigenstate below 200 GeV, which can be produced at the LHC. 
The prediction for the absolute Majorana neutrino mass scale m > 0.3 eV ensures 
that the model will be tested by future cosmological tests and (3(3q v searches. Rates for 
lepton flavour violating processes lj — > U + 7 typically lie in the range of sensitivity of 
coming experiments, with BR(/i — > erf) ^ 10 -15 and BR(r — > fij) > 10~ 9 . 

Finally, we mention that there have been attempts to realize the Harrison-Perkins- 
Scott mixing pattern in Eq. (151)1 at some high energy scale, and to correct its predic- 
tions by renormalization group evolution [82,83]. For a survey of related attempts see 
Ref. [84]. 

5.3 Absolute scale of neutrino mass and Offlfl 

Neutrino oscillations are blind to whether neutrinos are Dirac or Majorana. As we 
have seen, on general grounds, neutrinos are expected to be Majorana [9]. Neutrino- 
less double beta decay and other lepton number violation processes, such as neutrino 
transition electromagnetic moments [21,22] [85,86] are able to probe the basic nature 
of neutrinos. 

The significance of neutrinoless double beta decay stems from the fact that, in a 
gauge theory, irrespective of the mechanism that induces Qv(3(3 , it necessarily implies 
a Majorana neutrino mass [17], as illustrated in Fig. |§1 Thus it is a basic issue. 
Quantitative implications of the "black-box" argument are model-dependent, but the 
theorem itself holds in any "natural" gauge theory. 

Conventional neutrino oscillations are also insensitive to the absolute scale of neu- 
trino masses [18-20]. Although the latter will be tested directly in high sensitivity 
tritium beta decay studies [87], as well as by its effect on the cosmic microwave back- 
ground and the large scale structure of the Universe [88-90] QvfiP may give valuable 
complementary information. For example, as seen above, the model [79] gives a 
lower bound on the absolute Majorana neutrino mass m v ^0.3 eV and may therefore 
be tested in Ovfifi searches. 

2 There have been realizations of the A4 symmetry that also predict the solar angle, e. g. Ref. [80]. 
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Figure 9: Neutrinoless double beta decay and Majorana mass are equivalent [17]. 



Now that oscillations are experimentally confirmed we know that 0v(3(3 must be in- 
duced by the exchange of light Majorana neutrinos, the so-called "mass-mechanism". 
The corresponding amplitude is sensitive [19,20] both to the absolute scale of neutrino 
mass, as well as to Majorana phases [9], neither of which can be probed in oscilla- 
tions [18,19]. 

Fig. El shows the estimated average mass parameter characterizing the neutrino 
exchange contribution to Oz//?/? versus the lightest and heaviest neutrino masses. The 
calculation takes into account the current neutrino oscillation parameters in [14] and 
state-of-the-art nuclear matrix elements [91]. The upper (lower) panel corresponds to 
the cases of normal (inverted) neutrino mass spectra. In these plots the "diagonals" 
correspond to the case of quasi-degenerate neutrinos [79] [92] [93] In the normal hi- 
erarchy case there is in general no lower bound on the l3(3o v rate since there can be 
a destructive interference amongst the neutrino amplitudes. In contrast, the inverted 
neutrino mass hierarchy implies a "lower" bound for the f3f3o u amplitude. A specific 
normal hierarchy model for which a lower bound on /3/3o^ can be placed has been given 
in Ref. [80]. An interesting feature is that such lower bound depends, as expected, on 
the value of the Majorana violating phase 0i, as indicated in Fig. ^2 

The best current limit on (m u ) comes from the Heidelberg-Moscow experiment. 
There is also a claim made in Ref. [94] (see also Ref. [95]) which will be important to 
confirm or refute in future experiments. GERDA will provide an independent check 
of this claim [96]. SuperNEMO, CUORE, EXO, MAJORANA and possibly other 
experiments will further extend the sensitivity of current j3j3 0u searches [97]. 

5.4 Other phenomena 

If neutrino masses arise a la seesaw the dynamics responsible for generating the small 
neutrino masses seems most likely untestable. In other words, beyond neutrino masses 
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Figure 10: OvfiP amplitude versus current oscillation data, from Ref. [91]. 

and oscillations, theories can not be probed phenomenologically at low energies, due 
to the large scale involved. However when the breaking of lepton number symmetry 
is spontaneous there is a dynamical "tracer" of the mass-generation mechanism which 
might be probed experimentally. 

5.4.1 Majoron physics 

If neutrino masses follow from spontaneous violation of global lepton number, the exis- 
tence of the Goldstone boson brings new interactions for neutrinos and Higgs boson(s) 
which may lead to new phenomena. While neither of these is expected within the usual 
high-scale seesaw schemes, both could lead to detectable signals in low-scale models of 
neutrino mass. 

As already mentioned, the majoron may couple substantially to the SM Higgs boson, 
which can therefore have a sizeable decay branching ratio [36-39] into the channel in 
eq. ()39jl . Such "invisible" channel is experimentally detectable as missing energy or 



28 




\) 0.2 0.4 0.6 0.8 1 

I COS^j) I 



Figure 11: Lower bound on | (m ee ) |/Am| TM vs |cos(0i)| where <pi is a Majorana phase 
in the model of Ref . [80] . The lines in dark (red) and grey (green) correspond to normal 
and inverse hierarchy. 

transverse momentum associated to the Higgs [40,41]. Therefore in low-scale models 
of neutrino mass the neutrino mass-giving mechanism may have a strong impact in the 
electroweak sector. 

Majoron-emitting neutrino decays can also be lead to detectable signals in low-scale 
models of neutrino mass, 

f 3 ->• V + J, 

For example, if the neutrinos decay in high density media, like supernovae, charac- 
terized by huge matter densities, then the "matter-assisted" decays would lead to 
detectable signals in underground water Cerenkov experiments [98]. 

5.4.2 New gauge boson 

Although in the usual large-scale seesaw with gauged B-L symmetry, there are new 
gauge bosons associated to the neutrino mass generation these are too heavy to give 
any detectable effect. Only when the B-L scale is low, as in the model discussed in 
Sec. V>. 1.51 or the model considered in Ref. [99], there will exist a light new neutral gauge 
boson, Z' that could be detected in searches for Drell-Yan processes at the LHC. 

5.4.3 Lepton flavour violation 

In the presence of supersymmetry, seesaw phenomenology is richer. A generic feature of 
supersymmetric seesaw models is the existence of processes with lepton flavour violation 
such as fi~ — > e~7. Supersymmetry contributes through the exchange of charginos 
(neutralinos) and sneutrinos (charged sleptons) [100-102], as illustrated in Fig. IT2l 
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Figure 12: Supersymmetric Feynman diagrams for l { — > lj 7. 



Similarly the nuclear /1 — e conversion arises, as indicated in Fig. El The rates for 




Figure 13: Contributions to the nuclear u. — e conversion: (a) long-distance and (b) 
short-distance. For numerical results see Ref. [108] 

these process can both be sizeable. As an example, consider the rates for the a~ — > e~7 

decay, given in Fig. H3I 

The calculation leading to Fig. [TJ] is done in the framework of the supersymmetric 

double/inverse seesaw model in Sees. I3.1.11 and l3~2.3l [103], This allows one to analyse 

the interplay of neutral heavy lepton [104] 3 and supersymmetric contributions [102]. 

The figure shows the Br{ji — > ej) contours in the (M, -^)-plane (logarithmic scales) 

for hierarchical light neutrinos with mi = eV (left panel) and for degenerate light 

neutrinos with mi = 0.3 eV (right panel). The shaded contours correspond to (from left 

to right): Br(fi — > ej) = I0~ 15 '~ 13 '~ n '~ 9 . For large M the estimates are similar to those 

of the standard supersymmetric seesaw. However, if the neutral heavy leptons are in the 

TeV range (a situation not realizable in the minimal seesaw mechanism), the Br{^i — > 

cy) rate can be enhanced even in the absence of supersymmetry. This is indicated by 

the contours in the lower left, which depict the contribution from neutral heavy leptons 

only. For such M values around TeV or so, the quasi-Dirac neutral heavy leptons may 

3 Since in this model flavor and CP violation can occur in the massless neutrino limit, the allowed 
rates are unsuppressed by the smallness of neutrino masses [52,53, 104-106]. 
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Figure 14: Br(fi — > ej) in the supersymmetric inverse seesaw model of Ref. [103]. 

be directly produced at accelerators [107]. Note that in order to have low M then fi 
should be rather small, to keep neutrinos light. This is indicated by the diagonal lines, 
which indicate contours (top to bottom) of constant ji — 1, 10 -3 , 10 -6 , 10~ 9 GeV. The 
vertical lines are contours of Br(fi — > cy) in the standard supersymmetric seesaw. 

Similarly, the rates for the nuclear /i~ — e~ conversion in Fig. ED [108] fall within 
the sensitivity of future experiments such as PRISM [109]. 

5.4.4 Reconstructing neutrino mixing at accelerators 

Low-scale models of neutrino mass, considered in Sec. 13.21 offer the tantalizing pos- 
sibility of reconstructing neutrino mixing at high energy accelerators, like the "Large 
Hadron Collider" (LHC) and the "International Linear Collider" (ILC). 

A remarkable example is provided by the models where supersymmetry is the origin 
of neutrino mass [45], considered in Sec. 13.2.21 A general feature of these models is that, 
unprotected by any symmetry, the lightest supersymmetric particle (LSP) is unstable. 
In order to reproduce the masses indicated by current neutrino oscillation data, the 
LSP is expected to decay inside the detector [50] [110]. 

More strikingly, LSP decay properties correlate with neutrino mixing angles. For 
example, if the LSP is the lightest neutralino, it should have the same decay rate into 
muons and taus, since the observed atmospheric angle is close to 7r/4 [111-113]. Similar 
correlations hold irrespective of which supersymmetric particle is the LSP [114] and 
constitute a smoking gun signature of this proposal that will be tested at upcoming 
accelerators. 

There are other examples of low-scale models that nicely illustrate the possibility 
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tan 2 9 Atm tan 2 8 At-, 



Figure 15: LSP decays trace the atmospheric mixing angle [111], 
of probing neutrino properties at accelerators [115]. 



5.5 Thermal leptogenesis 

Now we briefly discuss one of the cosmological implications of neutrino masses and 
mixing, in the context of seesaw schemes. It has long been noted [116] that seesaw mod- 
els open an attractive possibility of accounting for the observed cosmological matter- 
antimatter asymmetry in the Universe through the leptogenesis mechanism [117]. In 
this picture the decays of the heavy "right-handed" neutrinos present in the seesaw 
play a crucial role. These take place through diagrams in Fig. in order to induce 
successful leptogenesis the decay must happen before the electroweak phase transi- 
tion [118] and must also take place out-of-equilibrium, i. e. the decay rate must be 
less than the Hubble expansion rate at that epoch. Another crucial ingredient is CP 
violation in the lepton sector. The lepton (or B-L) asymmetry thus produced then gets 




Figure 16: Diagrams contributing to leptogenesis. 

converted, through sphaleron processes, into the observed baryon asymmetry. 

In the framework of a supersymmetric seesaw scheme the high temperature needed 
for leptogenesis leads to an overproduction of gravitinos, which destroys the standard 
predictions of Big Bang Nucleosynthesis (BBN). In minimal supergravity models, with 
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m 3 /2 ~ 100 GeV to 10 TeV gravitinos are not stable, decaying during or after BBN. 
Their rate of production can be so large that subsequent gravitino decays completely 
change the standard BBN scenario. To prevent such "gravitino crisis" one requires an 
upper bound on the reheating temperature Tr after inflation, since the abundance of 
gravitinos is proportional to Tr. A recent detailed analysis derived a stringent upper 
bound Tr ^ 10 6 GeV when the gravitino decay has hadronic modes [119]. 

This upper bound is in conflict with the temperature required for leptogenesis, 
Tr > 2x 10 9 GeV [120]. Therefore, thermal leptogenesis seems difficult to reconcile with 
low energy supersymmetry if gravitino masses lie in the range suggested by the simplest 
minimal supergravity models. Their required mass is typically too large in order for 
them to be produced after inflation, implying that the minimal type I supersymmetric 
seesaw schemes may be in trouble. Two recent suggestions have been made to cure 
this inconsistency. 

One proposal [121] was to add a small R-parity violating XiV c iH u Hd term in the 
superpotential, where v c i are right-handed neutrino supermultiplets. One can show 
that in the presence of this term, the produced lepton-antilepton asymmetry can be 
enhanced. An alternative suggestion [122] was made in the context of the extended 
supersymmetric seesaw scheme considered in Sec. 13.1.51 It was shown in this case that 
leptogenesis can occur at the TeV scale through the decay of a new singlet, thereby 
avoiding the gravitino crisis. Washout of the asymmetry is effectively suppressed by 
the absence of direct couplings of the singlet to leptons. 

6 Non-standard neutrino interactions 

Most neutrino mass generation mechanisms imply the existence of dimension-6 sub- 
weak strength eGp non-standard neutrino interaction (NSI) operators, as illustrated in 
Fig- El These NSI can be of two types: flavour- changing (FC) and non- universal (NU). 
They are conceptually interesting for neutrino propagation since their presence leads 
to the possibility of resonant neutrino conversions even in the absence of masses [123]. 

NSI may arise from the non-trivial structure of charged and neutral current weak 
interactions characterizing seesaw-type schemes [9]. While their expected magnitude is 
rather model-dependent, it may well fall within the range that will be tested in future 
precision studies [68]. For example, in inverse seesaw model of Sec. 13.2.31 the non- 
unitary piece of the lepton mixing matrix can be sizeable and hence the induced non- 
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Figure 17: Flavour-changing effective operator for non-standard neutrino interaction. 

standard interactions may be phenomeno logically important. With neutrino physics 
entering the precision age it becomes a challenge to scrutinize the validity of the unitary 
approximation to the lepton mixing matrix, given its theoretical fragility [9]. 

Relatively sizable NSI strengths may also be induced in supersymmetric unified 
models [100] and models with radiatively induced neutrino masses, discussed in Sec. 13.2.11 

6.1 Atmospheric neutrinos 

The Hamiltonian describing atmospheric neutrino propagation in the presence of NSI 
has, in addition to the usual oscillation part, another term i^NSb 



Here +(— ) holds for neutrinos (anti- neutrinos) and e and e' parameterize the NSI: 
\^2GpNf€ is the forward scattering amplitude for the FC process + / — ► v T + f 



scattering (NU). Nf is the number density of the fermion / along the neutrino path. 

It has been shown [124] that in such 2-neutrino approximation, the determination of 
atmospheric neutrino parameters Am^ TM and sin 2 9 ATM is hardly affected by the presence 
of NSI on down-type quarks (/ = d). Future neutrino factories will substantially 
improve this bound [125]. 

6.2 Solar neutrinos 

Are solar oscillations robust? Do we understand the Sun, neutrino propagation and 
neutrino interactions well enough to trust current oscillation parameter determinations? 
Reactors have played a crucial role in identifying oscillations as "the" solution to the 
solar neutrino problem [62]. 




(52) 



and \f2GpNfS 1 represents the difference between + / and v T + f elastic forward 
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Thanks to KamLAND we have ruled out an otherwise excellent solution of the solar 
neutrino problem based on spin-flavour precession due to convective zone magnetic 
fields [126,127]. The absence of solar anti-neutrinos in KamLAND [128] has been used 
to establish robustness of the LMA solution with respect to spin-flavour precession due 
to convective zone magnetic fields. 

Thanks to KamLAND one could also establish robustness of the LMA solution with 
respect to small density fluctuations in the solar interior [129,130], as could arise, say, 
from solar radiative zone magnetic fields [131]. 

Finally, again thanks to KamLAND we have ruled out an otherwise excellent solu- 
tion of the solar neutrino problem based on non-standard neutrino interactions [132]. 
However, in contrast to the atmospheric case, non-standard physics may still affect 
neutrino propagation properties and detection cross sections in ways that can affect 
current determinations [133]. This implies that the oscillation interpretation of solar 
neutrino data is still "fragile" with respect to the presence of non-standard interactions 
in the e — r sector, though the required NSI strength for non-robustness to set in is 
quite large. In contrast, one can show that even a small residual non-standard inter- 
action of neutrinos in this channel can have dramatic consequences for the sensitivity 
to #i3 at a neutrino factory [134]. It is therefore important to improve the sensitivities 
on NSI, another window of opportunity for neutrino physics in the precision age. 
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